Key species and properties for perturbations of food webs 
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It is a challenge to predict the response of a large, complex system to a perturbation. Recent 
attempts to predict the behavior of food webs have revealed that the more complex the system, the 
more precisely the elements of the system must be measured. As a result, the amount of effort needed 
to understand a system grows quickly with its complexity. Here, we show that not all elements must 
be measured equally well, suggesting a more efficient allocation of effort to understanding complex 
systems is possible. We then develop an iterative technique to efficiently arrive at this solution. 
Finally, in our assessment of model food webs, we find that it is most important to precisely measure 
the mortality and predation rates of large, generalist, top predators. Prioritizing the study of such 
species will make it easier to understand the response of complex food webs to perturbations. 



I. INTRODUCTION 

Predicting the result of environmental perturbations, 
such as the arrival of new species or habitat change, is 
a major goal in ecology [l|, |38|. What makes this chal- 
lenging is the complex interconnected nature of ecolog- 
ical systems. In any densely-connected system, a per- 
turbation of one element can percolate across the net- 
work of interactions. This is particularly true for the 
complex food webs that form the backbones of most 
ecosystems [l^, [33|, [s^, [s^- Even perturbations acting 
only on a small subset of species may thus propagate 
through the network and lead to serious systemic changes 
[4, 9, 27, 31, 36, 45]. 

A central factor determining the response of a food 
web to perturbations is its topology, the precise map of 
predator-prey interactions. It has been shown that topo- 
logical properties affect local and global dynamical stabil- 
ity [il, M-, [HI, I28| - [3Q| and other notions of robustness such 
as the likelihood of secondary extinctions [7J. Often, the 
food web topology can indicate the relative importance 
of species for different applications [l, [l9|, |20|, [HI • How- 
ever, topology alone is generally not sufficient for reliable 
predictions [32]. For instance, measures taking into ac- 
count the biomass flows between species correlate much 
better with the results of numerical simulations than do 
topological indices alone [21]. 

In addition to knowledge of the topology, an accurate 
prediction of the impact of perturbations requires infor- 
mation about underlying biomass flows and the control 
coefficients characterizing the nonlinearity of processes. 
Such parameters require extensive measurements and er- 
rors in their estimation quickly reduce the accuracy of 
predictions about how food webs respond to perturba- 
tions [33|,|43. 

A lack of precise information on biomass flows and con- 
trol coefficients limits our ability to make precise predic- 
tions on the ultimate effect of perturbations. For in- 
stance, predictions for systems of more than 25 species 



are practically impossible, unless very detailed informa- 
tion is available ^ . 



Here, we ask if food web responses to a perturbation 
are more sensitive to particular species or parameters. If 
we can identify such influential elements in advance, we 
should be able to make more precise predications about 
dynamics with a given effort in paramaterization. 



In this paper, we investigate the predictability of re- 
sponses to perturbations in a broad class of food web 
models. Our results show that the different parameters 
of different species do not need to be measured with the 
same accuracy. We use analytical calculations and nu- 
merical demonstrations to show that it is possible to as- 
sign to each species a value that indicates the importance 
of precise knowledge about this species for the quality of 
the prediction. Furthermore, we demonstrate that this 
importance can be estimated reasonably well from im- 
precise information. Finally, we identify which of the 
parameters of these important species are crucial to the 
prediction of responses to perturbations. 



The paper is structured as follows: We start in Sec- 
tion [IT] by introducing a method for predicting the im- 
pact of given perturbations in a broad class of food web 
models. The method is illustrated in Section [iTT] with 
two examples. In Section HVl we then derive measures 
for species' influence on others and for their sensitivity 
to perturbations. In Section |Vl we test these predictions 
in a series of numerical experiments. The numerical re- 
sults illustrate a feasible strategy for field studies, where 
mathematical analysis and experimental measurements 
are used to iteratively improve predictions. In Section |Vll 
we use computer experiments and statistical association 
to determine which parameters and types of species are 
most important to measure. 
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II. IMPACT EVALUATION 

Consider a biological system described by a set of state 
variables Xi^ . . . , Xn denoting, for instance, the abun- 
dances of established species in a food web. The system 
is now subject to a perturbation that is characterized by 
another set of variables Yi , . . . , Ym , for instance denoting 
the abundances of newly arriving species. 

We assume that, in the absence of the perturbation, 
the variables Xi, . . . , Xat are governed by a set of ordi- 
nary differential equations of the form 



dt 



Xi = Ai{Xi^ . . . , Xat, Yi, . . . , Im), 



(1) 



where Ai is a function representing the right-hand- side of 
the differential equations. For instance, the generalized 
model for food webs ^Ji] which we use below, describes 
the dynamics of the populations Xi , . . . , Xjy by N dif- 
ferential equation of the form 

= Gi{X) + Si{X,) - Li{X) - Mi{X,), (2) 

where G^, L^, M^, and Si are unspecified functions de- 
scribing, respectively, the gain by predation (Gi), the 
loss by predation(L^), the loss due to natural mortality 
(M^i), and the gain by primary production (Si) of the 
focal species. 

Following Ref. 35, we consider the case where the un- 
perturbed system resides in a stable equilibrium X* and 
where the perturbation is characterized by a small and 
constant F*, such as a new species persisting at a low 
constant abundance in the ecosystem due to initially pos- 
itive growth or constant influx. 

Because the stationary abundance, X*, of a given es- 
tablished species i is dependent on the new species F*, 
we can regard it as a function X* = X*(F*). We can 
then define the impact lij of a perturbation variable Y* 
on a resident species abundance X* as the change of X* 
per unit 1^*, i.e. 



dx: 



oy: 



(3) 



where we used |o to indicate that the derivative is eval- 
uated in the limit of vanishing densities of the arriving 
species Y^*. In other words, the entries of the impact ma- 
trix lij state the loss of units of the established species 
i per unit of arriving species j that enters the system. 

In simple models, the impact can be computed by first 
solving ([T]) for the stationary solution X*{Y*) and sub- 
sequently computing the derivative in (j3]). However, for 
more than three species, the analytical computation of 
the stationary solution becomes prohibitively difficult. 
Furthermore, we seek a general solution, that does not 
depend on the functional forms in the model. For these 
reasons, the explicit computation of the stationary solu- 
tion is not possible. 



Computing the stationary solution can be avoided by 
recognizing that the stationary density of a resident 
species X* can be considered as an implicit function that 
is defined as the solution of the stationarity condition 
= 74^(X*,y*). Using a corollary to the implicit func- 
tion theorem [22|, we can then write the impact matrix 
as 



1 = -J-^K, 



(4) 



where the superscript —1 indicates the matrix inverse. 
The matrix J is the so-called Jacobian, which is defined 
as the derivatives of Ai with respect to the abundances 
of established species, i.e. Jij = dAi/dXj\^. And, the 
matrix K is a x M defined by Kij = dAi/dYj\^, 
where |* indicates that the derivative is evaluated in the 
equilibrium under consideration. The matrix K thus cap- 
tures the direct impact of an arriving species on an es- 
tablished species, which is quantified by the reduction in 
production (or respectively increase in mortality) of the 
established species per unit of the arriving species. To 
establish K prior to the arrival therefore requires infor- 
mation about the species the arriving species is likely to 
interact with. 

In summary, (|4|) establishes a relationship between the 
direct proximal impact of the arriving species K, the in- 
direct ultimate impact I and internal interactions among 
the established species J. Beyond the example of per- 
turbations caused by an arriving species, (j4|) applies to 
press perturbations on an established community in gen- 
eral. As a note of caution, we remark that this relation- 
ship holds up to linear order. The impact- approximation 
therefore remains valid only as long as the perturbation 
caused is reasonably small. 

Explicitly, ([4]) means that we can obtain the ultimate 
impact matrix as a function of the Jacobian of the es- 
tablished community and the matrix of the direct per- 
turbations. For instance, consider a perturbation due to 
a single new species of abundance Y{ arriving in a sys- 
tem of established species. Because there is only a single 
arriving species, the matrix K is a vector. For each estab- 
lished species the entry Ki denotes the changes to its 
biomass intake and loss rates that are directly caused by 
the new species, e.g. by feeding or competition. The im- 
pact vector / then describes the response of each species 
to this change induced by K as the perturbation propa- 
gates through the network. 

In the equations above, we refer to the steady state 
of the system, which seems to imply that information 
about this state is required. However, relationship (|4]) 
remains valid independently of the specific steady state 
under consideration. When the matrices are evaluated, 
the steady state appears only in the Jacobian which con- 
tains elements of the form Jij = d/dXj{dXi/dt)\^. For 
instance, in the generalized food web model, this leads to 
expressions such as dGi / dXi \ ^ . Because we cannot eval- 
uate this expression without further assumptions, it is an 
unknown quantity. However, we note that for any spe- 
cific system the expression is simply a number. In other 
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words, this means that the unknown derivatives appear- 
ing in the Jacobian constitute unknown parameters of 
the model. 

So far, we have recognized that the unknown deriva- 
tives can be formally treated as unknown parameters of 
the model. However, as such, these parameters are hard 
to interpret and are thus not suitable for an ecological 
discussion of the results. We solve this problem by using 
a slightly different parametrization, which is obtained ei- 
ther by a special normalization procedure [l3| or directly 
by the identity 



Name Interpretation 



Value(s) 



dGi 



G* d\og Gi 



X* dlogXi 



(5) 



which is true for G*, X* > (a condition that is generally 
met by definition; see Ref. for the special case of X* = 
0). 

The expression on the right-hand-side of ([5]) is a prod- 
uct of two factors that have a direct interpretation in 
most apphcations. The first factor is a per-capita rate. 
Such rates have the dimension of inverse time and can 
be directly interpreted as characteristic turnover rates, 
in this case, as the per-capita growth rate of the mem- 
bers of species i by predation on other species. 

The second factor in ([5]) is a logarithmic derivative. 
Such derivatives are also called elasticities and have been 
proposed originally in economic theory [l8| and sub- 
sequently in metabolic control theory jsf and ecology 
[iol . [43} . They can be estimated well from observa- 
tional data and interpreted straightforwardly. For ev- 
ery power-law, fix) = ^dx^, the logarithmic derivative 
is d\og f / dlogx = p, independently of A or x. Thus, 
for instance, any linear function has an elasticity of one 
regardless of the slope. For functions that are not power- 
laws the elasticity still provides an intuitive non-linear 
measure of the sensitivity in the steady state. For in- 
stance, for the well-known Holling type-II functional re- 
sponse the corresponding elasticity is 1 if evaluated in 
the linear regime and at saturation. 

In summary, the identity ([5]) allows one to break the 
partial derivative of the process in the steady state into 
two constant factors, describing the per-capita rate and 
the sensitivity of the process, respectively. These factors 
are, therefore, well-defined ecological parameters in their 
own right, which can be understood and discussed even if 
the steady state of the system is unknown. For food webs, 
this parametrization leads to the Jacobian matrix expres- 
sions given in [X{l3j . According to (j4]), the Jacobian that 
is thus parameterized can then be used to relate a per- 
turbation to its eventual impact. In the following, we 
use this approach to discuss the prediction of the impact 
of perturbations on food webs. We have so far focused 
specifically on the arrival of new species. Let us however 
remark that the applicability of (|4|) is not limited to this 
case, but remains valid for all perturbations that can be 
quantified by a set of variables F*. 

We note that the approach to impact taken in this 
paper is closely related to [35|. Our main methodolog- 



Scale parameters - defining the biomass fiows in the steady state. 

Oii Rate of biomass turnover in species i (0;1] 

Contribution of predation by i [0,1] 

to loss rate of species j 

Xi,j Contribution of species i [0,1] 

to the prey of species j 

Pi Fraction of growth in species i (producers) 

gained by predation 1 (consumers) 

Pi Fraction of growth in species i 1 — pi 

gained by primary production 

ai Fraction of mortality in species i (top pred.) 

resulting from predation (O5I] (others) 

di Fraction of mortality in species i 1 — ai 

not resulting from predation 
Elasticities - Sensitivities of interactions to state variables. 

7i Sensitivity of predation in species i [0.5,1.5] 

to i's prey density 

Xij Exponent of prey switching 1 (passive) 

Pi Exponent of closure in species i [1,2] 

(pi Sensitivity of primary production in species i [0,1] 

to the density of species i 

ipi Sensitivity of predation in species i [0.5,1.5] 

to the density of predators 

TABLE I: Generalized model parameters as defined in Ref. [T3I. 



ical contribution is to apply this approach to general- 
ized models. The advantage of generalized modeling is 
its high numerical efficiency, which enables a detailed 
and statistically sound numerical exploration. For the 
practical application to real world food webs, generalized 
models offer additional advantages. In contrast to half- 
maximum concentrations and maximal growth rates used 
in conventional models, all parameters of the generalized 
model are defined in the state observed in nature. The 
parameters can therefore be measured directly without 
requiring a fitting procedure. Furthermore, the parame- 
ters are defined in such a way that their estimation from 
noisy data converges maximally fast [l8|,[25|. In practice, 
this should lead to a higher accuracy. 

The formulation of the generalized model is straight- 
forward. Based on the results of [l3| even large models 
with tens of species can be set up in few hours. Using 
the equations derived here, the impact of different species 
can then be computed in seconds on a small laptop, us- 
ing a simple algorithm [3]. Once the model has been set 
up, integration of new data requires entering new numer- 
ical values. The computation of impact and importance 
therefore presents only a small additional effort to the 
field work needed to measure parameters. 



III. EXAMPLES 

For illustration, we consider two examples, a simple 
predator-prey system and the 10-species food web shown 
in Fig. [TJ While the predator-prey system is analytically 
tractable, the second example requires numerical calcu- 



4 































Sensitivity 


Influence 



Response 
10 



Sens. 
&Infl. 



0.1 
-0.1 



■10 



FIG. 1: Responses of a food web to perturbations of different 
species. Shown are species (circles) and predator-prey rela- 
tionships (arrows, in the direction of the biomass flow). In 
a-j), the species represented as a thick- lined circle is affected 
by a perturbation that increases this species' biomass intake 
and thus leads to an increase in its abundance. The color 
code denotes the response of the system to this perturba- 
tion. The perturbation of each species i is normalized to the 
species turnover rate ai, i.e. Ki = a^. In particular, the effect 
on the perturbed species is always positive, but the amount 
its abundance changes is determined by the surrounding sys- 
tem. In k) and 1) the grey-scale denotes the sensitivity and 
the influence of each species respectively, approximating its 
propensity to experience or cause large impact (details are 
given in Section llVi or see [C] for their definitions (|C3p and 



lations. 

As a first example, we consider a class of predator-prey 
models in which a predator of abundance Xi consumes 
a producer of abundance X2. A detailed treatment and 
discussion of the stability of this system in terms of the 
generalized model parameters can be found in Refs. [ll|, 

The Jacobian matrix of this system near the steady 
state is 



ai(V^-/i) ai7 



(6) 



where represents each species' turnover rate and a 
the relative loss of the producer due to predation (in- 
stead of natural mortality). Furthermore, <p denotes the 
elasticity (i.e. sensitivity) of primary production to the 
producer abundance, 7 the elasticity of predation to pri- 
mary producer abundance, tjj the elasticity of predation 
to predator abundance, and fi the elasticity of natural 
mortality to a species' own abundance. 

We now consider the impact of the arrival of a compet- 
ing predator in the estabhshed producer-predator system. 
It can be assumed that this new predator has a direct neg- 
ative effect on the primary producer but no direct effect 
on the estabhshed predator, such that the perturbation 



matrix K contains the entries Ki^i = and i^^2,i < 0. 
As shown in detail in [Bl the impact on the established 
predator is 



ai 
det J 



7^2, 



and the impact on the producer is 

J2 = ^(M-^)^2, 



detJ' 



(7) 



(8) 



where det J > represents the determinant of the Jaco- 
bian matrix. 

We see that generally the impact on the established 
predator is negative. This result is intuitive as the estab- 
lished predator is now in exploitative competition with 
the arriving predator. Of particular interest is the case 
where the established predator suffers from linear loss 
(/i = 1) and has an effect on the producer that scales 
linearly with the predator abundance (?/^ = 1, i.e. there 
is no interference between the predators). In this limit, 
det J approaches 0. For the impact on the producer, /2, 
the small denominator is compensated hy — /n^ which 
also vanishes [13] • However, the impact Ii on the estab- 
lished predator approaches infinity. This is a manifes- 
tation of the well-known competitive exclusion principle, 
which precludes the coexistence of the predators in this 
case [I2I, [13 • The infinite value indicates that the impact 
of the arriving predator is disproportionately strong, be- 
cause it leads to a loss of stability. Also in larger systems, 
such infinite values of the impact occur whenever the ar- 
rival of the new species causes the loss of local stability 
due to a bifurcation. 

The assessment of impact in a larger food web can be 
carried out analogously, but requires numerical compu- 
tations in which the generalized parameters are set to 
specific values. In the following, we consider generalized 
models described in detail in Ref. [l5|. The topo logy of a 
food web was generated using the niche model [4l5. For 
this topology we then drew the generalized parameters 
uniformly and independently from the ranges given in 
Tab. [H A specific example of a food web in this family of 
models is shown in Fig. [H 

The panels in Fig. [T] show the response of the food 
web to the perturbation of different focal species. In the 
context of a newly arriving species, these perturbations 
represent the situation where the new species has a direct 
effect only on the one focal species. Perturbations of some 
focal species have greater impacts on the food web than 
others. For instance, in Fig. [1^) most species respond to 
the perturbation of the focal species while in Fig. [1^), 
few species respond to the perturbation of a different 
focal species. 



IV. SENSITIVE AND INFLUENTIAL SPECIES 

Now we identify two types of species that are partic- 
ularly important for perturbations. Sensitive species are 
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FIG. 2: The average absolute impact on a species with given 
sensitivity if a species of given influence is perturbed. The 
impact increases strongly with the influence of the perturbed 
species and with the sensitivity of the species feeling the 
impact. Parameters: System size N — 50, connectance 
C = 0.04. 



easily perturbed by disturbances propagating through 
the web, and influential species have a strong effect on 
other species, when perturbed. In this section, we use 
intuitive reasoning to derive measures for the sensitivity 
and influence of each species. A more formal derivation 
is given in O 

Close to a steady state, the dynamical properties of a 
system are characterized by its dynamical modes which 
are characterized in terms of the eigenvectors and eigen- 
values of the Jacobian matrix. For a given matrix, there 
are generally two different sets o f eig envectors, which are 
called right and left eigenvectors . For each eigenvalue 
Xk of the matrix there is a corresponding right eigenvec- 
tor v^^"^ and left eigenvector w^^\ 

One can visualize dynamical modes as vibrations trav- 
eling through a drum when it is struck. Here the different 
modes correspond to different notes that are played on 
the drum. The right eigenvectors characterize the pat- 
tern of vibration when a specific note is played. Specifi- 
cally, the elements of the right eigenvector describe how 
strongly the respective area of the drum vibrates in that 
note. The same is true for the food web. In a stable 
steady state that is hit by a short (pulse) perturbation 
the right eigenvectors govern how the system returns to 
the steady state after the perturbation. 

Drummers know how to play different notes by strik- 
ing different parts of the drum. This is captured by left 
eigenvectors. Specifically, the elements of the left eigen- 
vector for a given dynamical mode describe how strongly 
the specific mode is excited when the drum is struck in 
a given area. Similarly, in the food web, the left eigen- 
vectors characterize the strength of a specific dynamical 
response when a given species is perturbed. 

Intuitively, one can think of each dynamical mode as a 
possible response of the system to a perturbation. More 
precisely, the right eigenvector denotes the impact of re- 
sponse (which species "feel the vibrations"), while the 
corresponding left eigenvector denotes the type of per- 



turbation that can trigger a particular response (which 
species needs to be perturbed to "play a given note"). 
For instance, consider the pair of a right eigenvector 
V = (1,2) and a left eigenvector w = (1,0). If a per- 
turbation affects only the first species, the direction of 
the system's response due this dynamical mode is in the 
direction (1,2); i.e., the second species changes twice as 
much as the first. 

The strength of a mode's response is determined by its 
excitability. The excitability is for a given mode given by 
the corresponding eigenvalue l/|A/c| of the inverse Jaco- 
bian matrix J~^, which has the same eigenvectors as J 
(see O for niore details). Intuitively, the eigenvalue Xk 
of a dynamical mode indicates a system's resistance to 
a particular perturbation. The impact of such a pertur- 
bation is, therefore, inversely proportional to this resis- 
tance. 

In the case presented here, we consider that a pertur- 
bation continuously excites the same dynamical modes 
(press perturbation). The impact is therefore the com- 
bined continuous excitation resulting from the perturba- 
tion of these dynamical modes. 

The potential impact that a species feels due to a given 
dynamical mode is the product of the mode's excitabil- 
ity and the component of the right eigenvector on this 
species. The potential impact from all modes is the 
sum over the contributions from the individual modes. 
We define the sensitivity of a species, denoted S'e, as 
the logarithm of this potential impact, where the loga- 
rithm brings the numerical values into a more manage- 
able range. For a more formal derivation, see O 

The potential impact that a species causes by exciting 
a given dynamical mode is the product of the mode's ex- 
citability and the component of the left eigenvector on 
this species. The potential impact caused through all 
nodes is the sum over the excitations of the individual 
modes. Analogously to the sensitivity, we define the in- 
fluence^ denoted /n, as the logarithm of this potential 
excitation. For a more formal derivation, see O 

For the example food web in Fig. [TJ the sensitivity 
and influence of each node in the food web are shown 
in panels k) and 1). Comparison with the impact for 
this food web confirms that sensitive nodes were often 
affected by the perturbations happening elsewhere and 
that direct perturbation of influential nodes had a strong 
effect on the network. 

The close relationship between the sensitivity, the in- 
fluence, and the impact is also confirmed in Fig. O On 
average, the impact of a perturbation on a species in- 
creases strongly with its sensitivity, and with the influ- 
ence of the directly perturbed species. 

Inspecting a series of sample food webs, we observe 
that the number of the very influential and very sensitve 
species in each food web is small. For instance, we find 
that on average for each web only 18% of all species lie 
in the upper 30% of the sensitivity or influence range for 
that web. 

In summary, knowledge of the Jacobian of a specific 
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FIG. 3: Prediction quality when species are measured suc- 
cessively with higher precision. Starting on the left, all food 
web species are known with low precision. Advancing to the 
right the measurment error is reduced for one species at a 
time until all nodes have been measured with high precision. 
The different curves refer to different strategies for selecting 
the next species for improving precision. For the dashed line 
(empty symbols), the next species is selected randomly. For 
the solid curves, we first evaluate the sensitivity or influence 
of each species based on current knowledge and then select 
the species with the highest value of either sensitivity or in- 
fluence. For the dotted curves we select the species with the 
highest product of influence and sensitivity (Sejui). Error- 
carrying parameters of each node: a, /i, 0, 7, a, /3, x- Other 
Parameters: Initial error of each parameter 10%, final error 
2%, connectance C — 0.04. 



food web enables us to predict the impact of specific per- 
turbations, and also allows to gain a more general un- 
derstanding of the species' sensitivity and influence with 
regard to perturbations of the network. The main chal- 
lenge for impact assessment is thus to collect the nec- 
essary data for constructing the system's Jacobian. We 
show in the following that high precision is required only 
for some species. 



V. ITERATIVE PARAMETER ESTIMATION 

Intuitively, to understand the dynamics of a system, 
one would need accurate measurements of the properties 
of the most influential and sensitive species. However, 
our notion of sensitivity and influence is itself derived 
from the Jacobian matrix. Therefore, we can't identify 
the species that we need information from without hav- 
ing information on those species. We propose two ways 
to address this dilemma. In this section, we propose 
an iterative strategy in which existing preliminary in- 
formation is used to estimate the impact and sensitivity 
of species. This preliminary assessment is then used to 
obtain improved parameter estimates on seemingly im- 
portant species. Once additional data on these species 
becomes available, they can be used to further improve 
the estimates of the impact and sensitivity of species, re- 



fining the process. Thus, a cycle is formed in which the 
necessary information for precise impact predictions is 
iteratively assembled. 

We explore the quality of impact prediction in a se- 
ries of numerical experiments. In each experiment, the 
task is to predict the impact of a random perturbation 
to a food web that was generated according to the pro- 
cedure described in Sec. IIIIl The food web generation 
determines values of the generalized parameters of the 
true Jacobian of the food web. In addition to this true 
Jacobian we generate an estimated Jacobian, in which we 
simulate measurement errors by drawing each generalized 
parameter from a log-normal distribution centered on the 
corresponding parameter value used in the true Jacobian. 
We then compute the true impact of the random pertur- 
bation, / based on the true Jacobian, and the estimated 
impact / based on the estimated Jacobian. The quality 
of the impact estimation is then evaluated as the cosine 
of the angle between the true and the estimated impact 
vector, which is computed as Q = / • //(||/|| ||/||), where 
• denotes the scalar product and ||/|| is the norm of the 
vector /. The computation yields values between Q = 1, 
indicating an exact match of estimated and true impact, 
and Q = — 1, indicating that the true impact for every 
species is the exact opposite of the prediction. 

We now introduce a numerical implementation of the 
iterative strategy described above. We consider numeri- 
cal experiments in which the knowledge of the Jacobian 
is initially poor, such that the generalized parameters are 
drawn from a lognormal distribution with a standard de- 
viation of 10% of the true value. We furthermore assume 
that additional empirical work can be carried out on spe- 
cific species that brings the error in all parameters of the 
respective species down to 2%. Our aim was to carry out 
the precise measurements in the order that leads to the 
most rapid increase in the quality of impact prediction. 

For the purpose of demonstration, we consider four 
different protocols: a) precise measurements are carried 
out in random order, b) species are measured in the order 
of decreasing influence, c) species are measured in the 
order of decreasing sensitivity, d) species are measured 
in the order of the decreasing product of sensitivity and 
influence. The choice of species to measure next, was 
always based on the estimated Jacobian that is available 
at the time. Thus, only information is used that would 
also be available in the real world at the respective time. 

The results shown in Fig. [3] demonstrate that estimat- 
ing influence and sensitivity of the species prior to each 
measurement strongly increases the accuracy of predic- 
tions. This is particularly pronounced if measurements 
focus on the species with the highest product of sensitiv- 
ity and influence. For instance, after measuring 20% of 
all species according to this protocol, we attain a qual- 
ity of prediction comparable to measuring 60% — 80% 
of all species when species are chosen randomly. Using 
the estimation of influence and sensitivity to focus ob- 
servational or experimental efforts can thus significantly 
reduce the amount of empirical work that is needed to 
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FIG. 4: Average quality of an impact prediction in the presence of measurement errors. The average quality of an impact 
prediction if one node with the specified sensitivity (left) or influence (right) in a food web is subject to a measurement 
error. The different datasets refer to errors in different parameters (c.f. Tab. [J). Other parameters are: system size N = 50, 
connect ance C — 0.04, and the standard deviation of the relative measurement error 10%. For higher C, the effect of x ^^nd f3 
increases (not shown). Noise around high importance and sensitivity values is due to the the relatively rare occurance of these 
values in the numerical experiments. 



achieve a given prediction quality. 



VI. MOST IMPORTANT PARAMETERS AND 
SPECIES TO MEASURE 

The iterative refinement procedure proposed above, 
needs some initial information on the system as a start- 
ing point. We therefore explore in this section, what 
types of parameters and what types of species were most 
important to measure. 

To get an initial intuition of the importance of different 
parameters for impact prediction, we consider a situation 
where the estimated Jacobian is identical to the true Ja- 
cobian except for a single parameter that carries an er- 
ror. Fig. [4] shows the quality of the estimated impact as 
a function of the influence and sensitivity of the species 
affected and the type of parameter varied. The figure 
furthermore shows the decrease in quality for sensitive 
and influential species depends on the parameter under 
consideration; precision in the elasticity of the mortal- 
ity /i, and of the elasticity of predation ijj with respect 
to predator abundance were the most important. Fig. [H 
thus confirms our intuition that not all parameters need 
to be measured to the same level of accuracy. 

To determine which species are most important to 
measure in the absence of knowledge about the Jacobian, 
we look for simple correlations between sensitivity or in- 
fluence and species properties (see below) in a set of 10^ 
stable food webs. Each sample web consists of TV = 50 
species and has a connectance of 0.04. They are sam- 
pled by generating a niche model topology, then drawing 
the generalized parameters from uniform distributions in 
their ranges given in Tab.H and retaining only stable con- 
figurations. For a given web, consider that each of the 
species i has a property Xi and sensitivity or influence 
Hi. Then we denoted the correlation coefficient between 



X and y diS R = {xi — x){yi — y)/ yx^^ where i runs 
over all species in a web, and where x denotes the mean 
value of this property. 

In the analysis, we consider the correlations (not causal 
effects) of sensitivity and influence with the following po- 
tential biological indicators: 

• Generality, or the number of prey species of the 
focal species. 

• Vulnerability, or the number of predators of the 
focal species. 

• A binary value that is 1 if the focal species is a 
primary producer and otherwise. 

• The trophic level TL, which we calculate by solving 
a set of linear equations, such that TLi = 1 for 
primary producers and TLi = l+mean(TLp^e2y) for 
consumers, where medin{TLprey) denotes the mean 
trophic level of z's prey. 

• The biomass turnover rate (generalized model pa- 
rameter a), indicating the amount of bimass an in- 
dividual consumes in comparison to its own mass. 

• The network degree of a species, i.e. the total num- 
ber of direct feeding interactions with other species. 

• The weighted topological importance of a species 
WP as introduced in detail in[20|. For each species, 
the value of WP indicates the indirect interac- 
tions from other species, based on the topology and 
biomass flows of a food web. The step parameter 
s indicates the maximum number of direct interac- 
tions, through which indirect effects are perceived. 
For instance, for s = 2, only effects on a species are 
considered from species that have a common prey 
or predator. 
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FIG. 5: Correlations of species' properties with their sensi- 
tivity and influence. The pattern of correlations is consistent 
with large predators playing an important role for a system's 
response to perturbations. 

The correlation analysis, presented in Fig. [5l shows 
that high trophic levels and low biomass turnover rates 
correlate strongly with sensitivity and influence. This re- 
sult suggests that top predators and other large species, 
despite their typically small total biomass, play a dispro- 
portionate role in the systems' response to perturbations. 
Without specific information on biomass flows, these are 
the species that should be targeted for initial parameter- 
ization. 

The sensitivity of a species is highly correlated with its 
generality, while its influence appeares to be independent. 
Intuitively, this might mean that a species with a large 
number of prey is likely affected if any of its prey species 
changes, while a perturbation of this species has little 
effect on each of the prey species. 

The weighted topological importance correlates 
strongly with sensitivity and weakly with influence. This 
might mean that indirect effects influence how perturba- 
tion affect a species. However, the indirect effects from 
distant species had little effect. 

In summary, high trophic levels and low biomass 
turnover rates correlate strongly with sensitivity and in- 
fluence, while generality and large indirect effects corre- 
late only with sensitivity. However, determining whether 
these correlations are causal requires statistical analyses 
beyond the scope of this paper. 



VII. CONCLUSIONS 



food webs. Our results show that there are typically 
a small number of species that are highly important, 
because they react sensitively to perturbations, have a 
strong influence on others, or both. 

While we have focused exclusively on food webs, we 
note that the same approach can likewise be applied to 
other networks of nonlinear interactions that are found in 
metabolism [39| , gene regulation [lo| , and cellular popu- 
lation dynamics |46j . 

Our results suggest that the potential impact of envi- 
ronmental perturbations on food webs can be predicted 
with reasonable accuracy if the most relevant parameters 
for only a small number of important species in the web 
are measured well. We have shown how to find these 
species iteratively, and how to pre-select potentially im- 
portant species based on their biological properties. In 
addition, we have identified the parameters that are par- 
ticularly important to measure accurately. 

For real world food webs, identifying important species 
requires some initial information about the system. How- 
ever, we demonstrated in numerical experiments that an 
iterative approach is feasible where preliminary inaccu- 
rate information is used to estimate the importance of 
different species, which can then be used to guide further 
field work, leading to refined predictions. Computing the 
importance of species based on a generalized model [l3| 
can be done in seconds, even for large systems with tens 
of species, by using a simple algorithm [3j. 

Our correlation analysis suggests that it is most im- 
portant to obtain precise parameter estimates for large, 
generalist consumers at top trophic levels. This analysis 
was not meant to determine which factors cause a species 
to be influential in a food web. Many of the biological 
traits we measured are highly correlated, and parsing out 
their independent effects was not a goal of our study. One 
could investigate hypotheses about biological traits that 
lead species to be more sensitive or influential in a food 
web, but that would require a different statistical ap- 
proach than used here and is best the subject of future 
work. 
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Previous work has suggested that without near per- 
fect information on a range of parameters, it may be in- 
tractable to predict the effects of perturbances to large, 
complex systems [35|. In this paper, we proposed a 
method to predict the impact of perturbations on com- 
plex systems more efficiently. We used this method to 
investigate the relative importance of different species in 



Appendix A: Generalized model Jacobian matrix 

In this appendix, we show the explicit expressions of 
the Jacobian matrix for general food webs of the form([T]) 
[l3|. The diagonal and off-diagonal entries of the Jaco- 
bian depending on the generalized parameters in Tab. [J 
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are seperately given by 



- {^3.i + X] /5/c,iA/e,j ((7/c - l)X/c,i)^ j . 

Intuitively, the terms in the first line of each expression 
result from biomass intake, either by primary production 
of the focal species z (p^ = = 0) or from feeding on 
a species j [pi = 0,pi = l,Xi,j 7^ 0)- The expressions 
in the second line of each expression are due to the loss 
of the focal species i either by natural mortality (aifii) 
or predation by species j {j3j^i ^ 0). The sums over k 
translate the indirect effects of apparent competition if i 
and j share a predator k. For a more detailed discussion, 
please refer to Ref. [l3|. 

We note that J characterizes the dynamics near any 
steady state. Stability of any particular steady state (cor- 
responding to a particular set of generalized parameters) 
is ensured by checking that all eigenvalues of the Jaco- 
bian have negative real parts [l6|, . 



Appendix B: Predator-prey example 

In this appendix we provide the explicit calculations 
of the simple producer-grazer example. In this system, 
a grazer of biomass density Xi consumes a primary pro- 
ducer of biomass density X2. The Jacobian matrix near 
the steady state is 



ai7i 



(Bl) 



A discussion of stability of this system in terms of 
the generalized modelling parameters can be found in 
Refs. [m, 

For calculating the impact of a given perturbation for 
each of the two species, we require the inverse of the 
Jacobian 

j-l ^ 1 I 0^2(^2 - Cr272 - ^2^2) -Q^l7l 

detJ I a20■2^pl (^li'^i - Pi) 



where 

det J = aia2(02 — Cr272 



^2P2){i^l - Pi) 



denotes the determinant of J. 

As a perturbation, we now consider a new grazer that 
has a direct negative effect on the primary producer, but 



no direct effect on the established grazer. The perturba- 
tion K is therefore given by i^i =0, K2 < 0. Now we 
plug this vector K in the impact equation (j4|) to obtain 
/ = — J~^(0, i^2)^- The impact on the grazer, /i, and the 
impact on the producer, I2 are therefore 



h = 



ai 
det J 



71 K2 and I2 = 



ai 
det J 



{pi-^i)K2. (B4) 



Appendix C: Derivation of influence and sensitivity 

Here, we derive the expressions for the influence and 
sensitivity of the food web species that are motivated in 
the paper with intuitive arguments. For this, we decom- 
pose the Jacobian matrix into the dynamical modes of 
the system. Then we use the impact-definition in (j4]) to 
identify sensitive species as those for which the expected 
impact of randomly excited modes is largest, and to iden- 
tify influential species as those for which the expected 
excitation of dynamical modes is largest when they are 
pertured. 

Each dynamical mode of the Jacobian matrix consists 
of a left eigenvector w''^\ a right eigenvector v'^^^ and a 
corresponding eigenvalue A/^. Writing all the right eigen- 
vectors as the columns of the matrix results in the trans- 
formation matrix V = {v^^"^ ... v^^"^). For the left eigen- 
vectors we analogously obtain the matrix W. In terms of 
these matrices, and of the diagonal matrix D containing 
the eigenvalues of J, we can write J = VDW-^. 

The inverse Jacobian matrix J~^, required for the im- 
pact calculation, has the same the same right and left 
eigenvectors v^^"^ and w^^"^ as the Jacobian matrix, corre- 
sponding to the eigenvalues l/Xk- Therefore can be 
written as J"^ = VD-^W^. 

Now we insert = VD~^W into the impact equa- 
tion dH). For a given perturbation vector K the impact 
is then 



w 



(k).K 



(Cl) 



where k runs over all eigenvalues and where • denotes the 
scalar product. Impact is thus the sum over the dynam- 
ical modes of the Jacobian matrix. 

We see that for a given given mode k in ()Cip the con- 
tribution to the impact on a given node depends on three 
factors. First, the product w^^^ • K determines which dy- 
namical modes are excited by the perturbation. Second, 
v^^^ determines which entries in / are affected by this dy- 
namical mode. And finally l/\k indicates how strongly 
the mode will be excited. 

For the sensitivity of a species, we do not want to ex- 
plicitly refer to any particular perturbation K. In ab- 
sence of additional information, we therefore consider the 
case of a perturbation affecting every dynamical mode k 
with identical probability. The sensitivity Sci can thus 
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be measured as the effects on node z, 




where denotes the absolute value of v^^"^ on node 

and where we used the logarithm to bring the numerical 
values into a more manageable range. 

For the influence of a species, we evaluate the impact 



of a perturbation affecting only this specific species K. In 
absence of additional information, the influence is mea- 
sured by the resulting excitation of all the dynamical 
modes 
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